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Abstract 

Classical 2D clock model is known to have a critical phase with 
Berezinskii-Kosterlitz-Thouless(BKT) transitions. These transitions 
have logarithmic corrections which make numerical analysis difficult. 
In order to resolve this difficulty, one of the authors has proposed the 
method called level spectroscopy, which is based on the conformal field 
theory. We extend this method to the multi-degenerated case. As an 
example, we study the classical 2D 6-clock model which can be mapped 
to the quantum self-dual ID 6-clock model. Additionally, we confirm 
that the self-dual point has a precise numerical agreement with the 
analytical result, and we argue the degeneracy of the excitation states 
at the self-dual point from the effective field theoretical point of view. 

1 Introduction 

Although low dimensional spin systems with a continuous rotational symme- 
try, like the 2D XY spin model, don't have a symmetry breaking phase at 
finite temperatures, it doesn't mean the absence of phase transitions. Ac- 
tually Berezinskii and Kosterlitz and Thouless 2| pointed out the vortex- 
antivortex excitation phase transition in the 2D XY model, which is known 
as the BKT transition. 
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The classical 2D clock model, which has a discrete rotational symmetry, is 
also known to have a BKT critical region. Jose et a/Q studied the 2D classical 
p-state clock model by the renormalization group analysis and pointed out 
the presence of the intermediate critical phase for large p. Elitzur et a/[4 3 have 
also obtained the same result through the correlation inequality argument. 

We map the 2D p-clock model onto its quantum ID p-clock model and 
analyze the latter. This is based on the fact that the (1 + d) dimensional clas- 
sical statistical mechanics can be mapped onto the d- dimensional quantum 
statistical mechanics. Here, we should be careful that the 2D classical clock 
model, introduced the following section, doesn't have the self-duality explic- 
itly, although its quantized Hamiltonian has the self-duality. It may seems 
strange, but thinking of the universality class; that is, the Z p Villain model 
and the 2D classical p-clock model are considered belong to the same univer- 
sality class, the self-duality of the quantum clock model is acceptable. Ad- 
ditionally, the corresponding continuum model(the dual sine-Gordon model) 
also has the self-duality. We will see an excellent numerical agreement with 
the self-duality. 

Generally, the BKT transition points were difficult to determine numer- 
ically because of the logarithmic corrections and the anomalous divergence 
of the correlation length. But the level spectroscopy method[5jjn] has been 
proposed to resolve these problems. Near a BKT transition point, some scal- 
ing dimensions change from relevant to irrelevant or vice versa. Choosing 
some applicable scaling dimensions, we can determine the BKT transition 
point by the crossing point of these scaling dimensions. In addition, we can 
eliminate logarithmic corrections accompanying the BKT transition. This is 
the brief explanation of the level spectroscopy. 

So far the level spectroscopy has been applied to the three cases. The first 
case is between the BKT critical phase and the phase of the non-degenerate 
ground state [S] The second is between the BKT critical phase and the 
phase of 2-fold degenerate ground state jS]. The third is between the BKT 
critical phase and the multi-degenerated ground state phase. Recently Tone- 
gawa et a/jS] studied the 1/3 plateau problem of the S = | antiferromagnetic 
XXZ chain with the Z 3 symmetry breaking and Otsuka et al[?} studied the 
2D AF 3-state Potts model with the Zq symmetry breaking case. 

In this paper, we apply the level spectroscopy to the ID quantum 6-clock 
model with duality and use the duality relation to check the level spec- 
troscopy results between the BKT critical phase and the multi-degenerated 
ground state from another point of view than former studies. In addition, we 
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discuss the degeneracy at the self-dual point from the field theoretical point 
of view. The self-dual point of this model is trivial but it will be important 
in the 2D AF 3-state Potts model which have no explicit self-duality. 



2 Theory 

2.1 Discrete model 

In this subsection, we review the ID quantum p-state clock model with a 
nearest neighbor interaction as a Hamiltonian limit of a 2D classical p-state 
clock model on a square lattice. The 2D classical Hamiltonian is 



n = -/?^{j,[cos(e M -e mi:; )-i] 

(i,j) 

+J T [cos(6 fJ - - 9y +1 )]}, 



(1) 



where Q^j = (2irr/p), r = 0, 1, . . . ,p — 1 is the clock spin variable. When we 
take the Hamiltonian limit (3J S — > 0, f3J T — > oo with A fixed [ID] . we obtain 
the ID quantum Hamiltonian |%]|1(J| 



N 

H = -2^2 |Acos(9 n - 6 n+ i) + cosp n | 



(2) 



n=l 



where A is the transverse field of the system. It can be interpreted as the 
inverse of the temperature. 

When we define the operators cr n ,T n as a n = expi9,r n = expip, we 
obtain another representation of the ID quantum Hamiltonian (J2J) as 



h = - i x K^+i + < ^n+i) + r n + r+} 



n=l 



where u = exp(2m/p), and 
/ 1 
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These operators satisfy (cr„) p = (T n ) p = 1. 

This model has a Z v charge symmetry. The Z v charge operator is defined 

as 

L 

u Q = U r„. (4) 

n=l 

which commutes with the Hamiltonian J3J). Uq has eigenvalues uj®,Q = 

0. 1, • • • ,p — 1. Because of the conservation, Uq splits H into the charge 
sectors whose corresponding eigenspaces are labeled by Q — 0, 1, • • • ,p — 1, 

1. e. H is able to be block-diagonalized. Also we can block-diagonalize it with 
the toroidal boundary conditions: 

a N+1 = exp [—Q \<Ti,Q = 0,l,...,p — l. (5) 

Q = is the periodic boundary condition and Q = p/2 (p even) is the twisted 
boundary condition. Hereafter we will assume that p is even. Implementing 
the twisted boundary condition to the system corresponds to adding a half- 
charge to the sine-Gordon model [7] JT] which is the continuum limit of the 
discrete model. 

There exits a unitary operator such that a n is transformed into and 
also T n is transformed into a n . For example we can take its matrix elements 

as 

I7y = 4=w (1 -° (W) . (6) 

This operator U gives the relation 

Ua n U + = r+ (7) 
UT n U + = a n . (8) 

By this relation, we get another representation of the ID quantum Hamilto- 
nian as 

L 

H ' = -Y,i X ( r « r -+l + r " r "+0 + ^ + O ■ ( 9 ) 

71=1 
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In this case, notice that the Z p charge operator should be 



(10) 



71=1 



The form Q, fllOj) is more useful for numerical calculation than the pribious 
form ©, when we choose the a n diagonal representation. 

There is another non-local transformation between a and T: 



On - 




— Tn+lT 


CTJV - 


-> &N 


— L N 


r„ - 


- fn 


n 






1=1 



(11) 



Applying this transformation to the Hamiltonian Q and the help of 

n n+1 



r r+ 

- 1 - n 1 - n+1 



ITdl 

i=l z=i 



(12) 



we get the dual form of (0) as 



£ [ A ( CT n + *») + r„r+ +1 + r+r rt+1 ] . 



(13) 



71=1 



Therefore the ID p-state quantum clock model satisfies the following duality 
relation, 



H(X) = XH(1/X). 



The self-dual point is A c 



The charge conjugation operator is defined as 

/ 1 



Cn i C-r, 



71=1 





1 
\0 1 



1 



0/ 



(14) 



(15) 
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which has eigenvalues ±1 and transforms a n and T n as 

c„a n c+ = <7+ (16) 

By this relation, it is easy to see that the Hamiltonian (jSJ) and (JHJ) commutes 
with C. 

Z p charge operator Uq satisfies 

CU Q C = U+. (18) 

So the eigenvalues uj® of Uq become u; - ^, that is, generally the eigenstates of 
Uq are not those of C. But the eigenstates of Q = 0,p/2 are invariant under 
this transformation, therefore these eigenstates have eigenvalues C = ±1. 
We will use this fact to classify the eigenstates. 



2.2 Effective model 

As an effective theory of a ID quantum spin system, the sine-Gordon model 
has been studied. Here, we will explain the Z p dual sine-Gordon model|T2]. 

where K is also inverse of temperature as A, a is an ultraviolet cutoff and (f) 
and 6 are periodic: 

<P = 0+^=, (20) 

e = e + ^. (2!) 

The field 6 is dual to cf> and they satisfy the relation 

d x <j)=-d y (K6), d y (f> = d x (Ke). (22) 

This model is the effective model for the Hamiltonian (jHj). (jUj). 

The vertex operators are defined as 

O m ,n = exp(imV26) exp(inV2<f)) . (23) 
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When = y e = 0, its scaling dimension is 

x m ,n{K) = \(n*K+?pj. (24) 

In a high temperature region (K small) the second term of (fTTH) is more 
relevant than the third term. So the renormalization group behaviour in the 
high temperature region can be studied through the Lagrangian 

£ hlgh = ^(V0) 2 + ^cosV20. (25) 

When K = 4 this sine-Gordon model shows the BKT transition, and that is 
studied by the level spectroscopy [7]. On the other hand, in the low temper- 
ature region the second term of (J19)) is more irrelevant than the third term. 
Then we deal with the sine-Gordon Lagrangian 

£io W = ^(W) 2 + 77^1 cospy/26 (26) 

27T 27TCT 

with the Z p symmetry breaking term. Here we got the first term through 
substituting ifTOjl to p2|) . 

Under the change of the cutoff a — > e a and parameterizing K/4 = 
l + |/o/2, the renormalization group equations for (f23|) are 

d V0(l) 2/ A d VS) /A /A / 0V N 

— ^— = -tfjxih —J— = -yA l m(i) (27) 

where I is related to e l — L for the finite system size. For the case (J26j) . by 
parameterizing K~ l p 2 /4 = 1 + y /2 and using yg instead of y^, we get the 
renormalization equations 

In the high temperature limit A — > (or K — > 0) this model has a non- 
degenerated ground state. The transition point between the non-degenerated 
ground state phase and the BKT critical phase can be determined by the 
K = 4 version of the level spectroscopy. In this case, the renormalized 
scaling dimensions are given in In the low temperature limit A — > 00, 
the model has a p-fold degenerated ground state. In this case, the simple 
K = 4 version of the level spectroscopy is no longer available. The BKT 
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Table 1: Operator content of the K = p 2 /4: sine-Gordon model with Z. p 
symmetry breaking term. TBC is short for twisted boundary condition. 
Others are under periodic boundary condition. 



field 


quantum number 


B.C. 


scaling 


renormalized 










dimension 


scaling 




Q 


c 






dimension 


exp(im^/26) 


m ytz p/2,p 




PBC 


Xm,0 




cos(V2p0/2) 


p/2 


1 


PBC 


x p/2,0 


\ (i - bo) 


sin( v / 2^/2) 


p/2 


-1 


PBC 


x p/2,0 


\ (i + iw) 


marginal 


P 


1 


PBC 


%p,0 









1 


PBC 




1 (i - bo) 


003(^0/2) 





1 


TBC 


u, 2 


1 (i - bo) 


sin(V2<£/2) 





-1 


TBC 


U >2 





critical transition point can be determined by the K = p 2 /4: version of the 
level spectroscopy (see also and [?]). The renormalized scaling dimensions 
of near K = p 2 /4 are shown in Tabl^note that implementing TBC is the 
same as considering the half-charge of the effective model [TT]). 

Assuming conformal invariance, the scaling dimensions x m>n are related 
to the energy gap of the finite size system with periodic boundary conditions 
as 

x m ,n = {E m<n {L) - E g (L)) , (29) 

where L is the system size, E g (L) is the ground state energy, and v is the 
velocity of the system [T5]. And the central charge c, which we will use 
to confirm the universality class of the system, is given by the finite size 
correction of the system as 

7TVC 

Eg{L) =e g L-— (30) 

where e g is the bulk ground state energy per site|14j. 

At last, we discuss the dual transformation in this dual sine-Gordon 
model, corresponding with the discrete case (Jll|) . (|14|) . When we apply the 
dual transformation 

(f)^p6, pB^cj) (31) 
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to the dual sine-Gordon model (|19p , we get the dual form 

2 

£duai = 7^77 (V#) 2 + -^cospV26 + -^- C osV2(f>. (32) 

When we substitute (J22J) to (jS2]) ) we obtain the same form as (fTTJj) aside 
the coefficients; assuming = ye and p = K, including the coefficients, we 
get the same equation as f)19|) . This is the self-dual point of the dual sine- 
Gordon model. Note that with the dual transformation ()31|). fields in Table 
^are exchanged as cos \^2p9/2 <-> cos \/2p/2 etc. These exchanges mean the 
degeneracy of the excitation spectra at the self-dual point. 

The quantum discrete model has the self-dual point (pHj ). and now this 
self-duality is extended to the continuum model. In the following section, we 
numerically estimate the self-dual point with good precision. 



3 Numerical Results 

In the previous section we have argued the general p-state clock model. In 
this section we show the numerical results for the p = 6 case. The ground 
states with the various Z p charges and the boundary conditions are obtained 
by the Lanczos method. The system sizes are up to L = 10. 

First, we determine the two BKT transition points. These points have es- 
sentially different property, because one is located between the non-degenerate 
ground state phase and the BKT critical phase, the another is located be- 
tween the 6-fold degenerate ground states phase and the BKT critical phase. 

Next, to confirm the universality class of the BKT critical phase, we 
calculate the ratio of the scaling dimensions, the central charge c and the 
parameter K of the sine-Gordon model. 



3.1 The BKT transition point in lower A 

The BKT transition point between the non-degenerate ground state phase 
and the BKT critical phase is given by the level crossing of the following 
energy gaps [I]: 

X2, = -^A£(Q = 2,PBC), (33) 
Znv 

x i = -^—AE(Q = 0,TBC,C = 1). (34) 
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0.351345 



0.351340 



L=10 with PBC 
L=10withTBC 



0.78231 



0.78232 



0.78233 



Figure 1: The crossing point for L = 10 in lower A 



AE is an energy gap between the ground state and the excitation state with 
the quantum number Q and a periodic boundary condition(PBC) or a twisted 
boundary condition(TBC). In terms of the toroidal boundary condition^, 
PBC corresponds to Q = 0, and TBC corresponds to Q = p/2. The ground 
state has Q = 0, C = 1 with PBC. 

Figure^ shows the crossing point for L = 10. On the BKT transition 
point, we can eliminate logarithmic corrections including higher terms by 
the level spectroscopy, but still we have finite size corrections due to the 
irrelevant field L^L^^ix = 4) ^3] which behaves order of 1/L 2 , 1/L 4 ,etc. 
In figureEl we show the extrapolation of the level crossing as L tends to 
infinity. The result is A = 0.78183. 



3.2 The BKT transition point in upper A 

The upper BKT critical point between the 6-fold degenerate ground state 
phase and the critical phase is given by the crossing point of the following 
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Figure 2: The lower BKT transition point. The extrapolated value is A 
0.78183. 
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Figure 3: The crossing point for L = 10 in upper A 



scaling dimensions [B] [?] 
9 

7 X 3,0 



9 L 



x 



A2nv 



AE(Q = 3,PBC,C= 1), 



0, TBC, C 



(35) 
(36) 



For Q = 0, 3, the charge conjugation C can be a good quantum number as 
mentioned at the end of the subsection 12.21 

The crossing point is shown in FigureI2]for L = 10. This crossing elim- 
inates the logarithmic correction order of 0(1/ In L). But still higher order 
corrections like 0(l/(lnL) 2 ) can remain in this case. Anyway we extrapolate 
BKT transition point for large limit of the system size L as in figureEJ We 
obtain A = 1.2851. 

The product of A U p per and Ai ower can be used to check the consistency of 



the self duality; that is, A upP er 
numerical data, we get 



• Ai OW er should be unity. Actually, from the 



A 



upper 



Ai ower 



1.0047. 



(37) 
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Figure 4: The upper BKT transition point. The extrapolated value is A = 
1.2851. 
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Figure 5: Level crossing occurs exactly at the self-dual point A = 1. These 
levels are AE(Q = 3, PBC, C = 1) and AE(Q = 0, TBC, C = 1) for L = 10. 
The other system sizes also have the exact crossing. 



3.3 Universality class 
3.3.1 Duality 

We numerically calculate the change of K via A by the ratio of scaling di- 
mensions: 

3,PBC,C 



AE(Q 



1 



2^3,0 



AE{Q = 0,TBC,C = 1) 



36/K 2 . 



(38) 



This ratio has the logarithmic correction 0(l/(lnL) 2 ) in the vicinity of the 
BKT transition points. But at the self-dual point, AE(Q = 3, PBC, C — 1) 
and AE(Q = 0, TBC, C — 1) exactly degenerate; therefore exactly K — 6. 
This is because the dual sine-Gordon model has higher symmetry at the self- 
dual point as was discussed in (p?Tj) . (JH^j) . We show this in figured The 
energy excitations: AE(Q = 3, PBC, C = 1) and AE(Q = 0, TBC, C = 1) 
exactly cross at A = 1. 



14 





& 1- 




3.9 - 




3.8 - 






W 




<\ 


3.7 - 




3.6 - 


w 


<\ 






3.5 - 




3.4 - 




3.3 - 




0.5 0.6 0.7 0.8 0.9 1 



A 

Figure 6: Ratio of scaling dimensions a^o/^i.o- m the critical region the 
ratio X2,o/^i,o — 4 



von Gehlen and Rittenbergjin] pointed out that the spectra of the Z p - 
symmetric quantum chains at the self-dual point are p-fold degenerated. Our 
novel point is not only the degeneracy but the connection of the discrete Z p 
model and the Z p dual sine-Gordon model. 

Even though when one study some other lattice model which doesn't 
have the self-duality explicitly, but if it could be mapped to the dual sine- 
Gordon model, one can check the self-dual point of the dual sine-Gordon 
model through this level crossing. For example, such an approach will be 
useful in the vicinity of the multicritical point studied by Otsuka et a/[?]. 



3.3.2 Ratio of scaling dimensions 

The ratio of scaling dimensions 

x 2 , AE(Q = 2,PBC) 



zi,o AE(Q = 1,PBC) 



4 + (L ) + higherorder 



(39) 
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is also useful to check the universality class of the BKT critical region. The 
finite correction behaves order of 0(L~ 2 ). Figure El shows the ratio. 

3.3.3 Central charge 

Generally, in a BKT critical region a renormalization equation flows to the 
Gaussian fixed line. Therefore the universality class of the BKT critical 
region agrees with the Gaussian model. The Gaussian model is known to 
have the central charge c = 1, so the BKT critical region is also characterized 
byc=l|inj. 

From the conformal field theory the central charge is related to the ground 
state of the finite system asjTT] 

W = ^-^( 1 + (ieW)) (40) 

where e g is a free energy per site and v is a velocity of the system. This is 
one of the important result of the conformal field theory. In this case the 
logarithmic correction is small enough, so the central charge is numerically 
a good index of universality class. 

We calculate the effective central charge for the system sizes up to L = 10 
and extrapolated to L — > oo. In Figure we show the extrapolated values 
of the central charge. This result is consistent with the two transition points 
we have determined. 



4 Conclusion 

We extended the level spectroscopy to be able to determine the BKT critical 
point between the multi-fold degenerated phase and the critical region. For a 
physical application we studied the ID quantum p-clock model which has the 
Z p symmetry. The numerical calculations were performed for p = 6. And also 
we discussed about the self-duality of the dual sine-Gordon model. In the sec- 
tion III-C, the self-dual point could be numerically determined without loga- 
rithmic corrections. This is because the level crossing of AE(Q = p/2, PBC) 
and AE(Q = 0,TBC) is the exact result from the duality (jlljl . There is no 
correction term in the langage of the dual sine-Gordon model (|19|) . On the 
other hand, the BKT transition point that we determined have some correc- 
tion beause we ignore irrelevant terms to derive (J25|) and (|26j) from (|19|) . Our 
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Figure 7: The central charge extrapolated from the finite system's results. 
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present approach will be useful for the other models which need more highly 
accurate calculations, for instance the cross-over near the multicritical point 
dealed in [?]. 
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